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Abstract 

Low energy scattering and bound state properties of the ttN, tttt and Ktt 
systems are studied as coupled channel problems using inversion potentials of 
phase shift data. In a first step we apply the potential model to explain recent 
measurements of pionic hydrogen shift and width. Secondly, predictions of 
the model for pionium lifetime and shift confirm a well known and widely 
used effective range expression. Thirdly, as extension of this confirmation, we 
predict an unexpected medium effect of the pionium lifetime which shortens 
by several orders of magnitude. The a meson shows a narrow resonance 
structure as a function of the medium modified mass with the implication of 
being essentially energy independent. Similarly, we see this medium resonance 
effect realized for the Ktt system. To support our findings we present also 
results for the p meson and the A(1232) resonance. 
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I. INTRODUCTION 



It is difficult to name hadronic systems for which the underlying QCD structure is man- 
ifest in low and medium energy observables ]IJ. The reason for this are the significantly dif- 
ferent masses of nucleons and mesons as compared to the quark masses. As a consequence we 
may observe the QCD effects at very short interhadronic distances with an energy scale for 
excitations in the GeV and not as it is customary for nuclear physics in the MeV region. In 
the language of quantum mechanical radial wave functions we anticipate a situation similar 
to atomic physics where it is difficult to disclose effects of the nuclear realm. The interplay 
between low energy nuclear physics and QCD is tough since the phenomena are still difficult 
to identify and to associate with a particular substructure. Triggered by the need of a scalar 
isoscalar medium weight a meson to understand the medium range nucleon-nucleon (NN) 
attraction, there exist many attempts to localize a resonance in L = 0, T = tttt scattering 
with a mass m a ~ 300 — 700 MeV Despite all attempts, experimental data do not 

support the anticipated narrow resonance and theory gets acquainted with this situation. 

We have been triggered by the development of a one solitary boson exchange model 
(OSBEP) |4[] for NN scattering to study the tttt and 7r-nucleon (ttN) system in terms of a 
potential model. The purpose of this investigation is to find a quantitative first guess for a 
dynamical model to describe the resonant particles p, a, A etc. It is intended to treat them 
explicitly as interacting multi pion and nucleon systems in an NN interaction model above 
meson production threshold. Undoubtedly, such a model is a subtle many-body problem 
for supercomputing. 

Despite this motivation, there exists considerable interest in studies of di-hadronic sys- 
tems which form Coulomb bound states [|5]J|. Precision measurements of shifts and widths 
of such hydrogen like two body systems are supposed to sense with high accuracy the strong 
interaction at low energy. The shift is mainly caused by a modification at short distances 
where both the Coulomb and hadronic interactions are present while the width results from 
decay into energetically open channels. The decay can lead to a lower lying state of the same 
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system or into reaction channels. Candidates for such studies are all oppositely charged 
hadron pairs. From this sample, we selected pionic hydrogen (A np ) as a bound 7r~-proton 
system and pionium (A nn ) as a bound it~tt + system both in relative s-states. 

To describe such atomic systems it requires to know the hadronic interaction a priori or 
make a fit with a potential ansatz. The hadronic interaction is known, in the sense that the 
partial wave phase shifts are sufficient to determine the interaction by quantum inversion. 
This mathematical method uses spectral theory to transform the boundary conditions in 
form of the S-matrix or Jost function into a local energy independent potential @. The 
radial Schrodinger and Klein-Gordon equations are equivalent to a Sturm-Liouville problem 
on the half axis which is central for Gelfand-Levitan-Marchenko inversion. In Sec. II we 
show the salient features of inversion and give in Sec. Ill the potential results for tcN, tttt 
and Kit scattering. With these potentials all ingredients are specified for the treatment 
of A np and A n7T as coupled channel problems and predictions are given in Sec. IV. These 
predictions are in good agreement with experiment and lent support to our potential model. 

As extension of the coupled channel potential model we treat 7riV, tctt and Kir scattering 
in the continuum. The coupled equations, with Coulomb effects included, confirm that 
isospin is hardly broken and that resonance properties of all systems are quantitatively well 
described. 

The p-wave resonances A and p are contained in Sees. III.B and D respectively. S-wave 
channels of and Kit scattering are treated in Sec. V. We attempt to solve the puzzle 
why a resonance is not visible in the phase shifts but at the same time meson exchange 
models require a medium weight a meson with definitive particle properties. The solution 
to this investigation is a medium modification of the two body scattering which is manifest 
only in the T = tttt channel and T = 1/2 Ktc channel. The medium modification is 
identified with an effective mass in the coupled equations leading to a resonance structure 
as a function of effective mass and not of energy. Thus we observe a a resonance width r w 1 
MeV practically independent of energy. Associated with this o resonance is a shortening of 
lifetime of pionium A nn by more than three orders of magnitude. The conclusions in Sec. 
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VI comprise the prediction of a very short range attraction close to the origin and a narrow 
high repulsive barrier at a relative distance r « 0.2 fm for the A, p, a and Kq resonances. 
Comparable short range repulsions are seen for other channels. A confirmation of the long 
range interaction in terms of meson exchange is given and most importantly a narrow a 
meson is predicted as a medium effect. 

II. POTENTIALS FROM INVERSION 

Contrary to the direct path to obtain a potential for elementary particle scattering from 
QCD or other microscopic models we apply quantum inversion to experimentally determined 
phase shift functions as input in Gelfand-Levitan-Marchenko equations |7j]. Nowadays, the 
inversion techniques for nucleon-nucleon scattering have evolved up to almost perfection 
for scattering data below pion production threshold. Numerically, input phase shifts can 
be reproduced for single and coupled channels with a precision of 1/100 of a degree, which 
is much lower than the experimental uncertainty. This accuracy and the possibility to 
test the inversion potential online, i. e. inserting the potential into the scattering equation 
and reproducing the input phase shifts, makes them a reliable and easy-to-handle tool 
for quantitative medium energy nuclear physics. Guided by this spirit, the utmost aim of 
quantum inversion is to provide a simple but accurate operator to reproduce data. This 
paradigm, however, proscribes to include sophisticated momentum dependencies or non- 
localities in the potential since this requires more information then can be extracted from 
phase shifts in a very limited energy domain. Our studies and applications are limited 
to the real phase shift domain but, for mathematical reasons, they are with smooth real 
functions extrapolated towards higher energies and infinity. This extrapolation is a kind 
of regularization which does not introduce spurious low energy phenomena and thus has 
no effect upon results and conclusions. This aspect of inversion has been investigated in 
important circumstances and we assume to be on save ground also in the applications studied 
hereafter. Furthermore, we are facing only single channel situations without Coulomb effects. 
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The basic equation of inversion is the Sturm-Liouville equation 

d 2 



dx 2 



+ q(x) 



y (x) = Xy(x). 



(1) 



We use the equivalent radial Schrodinger equation 

tpl{k,r) = k 2 tp e (k,r) 



d 2 £(£ + 1) 2/ic 2 . 



< r < oo 



(2) 



dr 2 r 2 (he) 2 

where Vg(r) is a local, k-independent operator in coordinate space and the factor 2/zc 2 / (he) 2 
guarantees the correct units. Boundary conditions for the physical solutions are 



and 



hmi^Ak, r) = 
i — >o 



£n 

lim ipe(k, r) = exp(i5e(k)) sin(A;r — — + 5e(k)) 



(3) 



(4) 



The Marchenko and the Gelfand-Levitan inversion are two inversion algorithm for the 
Sturm-Liouville equation which are briefly sketched for single channels and the case without 
a Coulomb reference potential. More details can be found elsewhere [0.f|. 

A. Marchenko Inversion 

The experimental information enters in the Marchenko inversion via the S'-matrix, which 
is related to the scattering phase shifts by the relation 



Si(k) = exp(2i5i(k)). 



We use a rational function interpolation and extrapolation of real data 5e(k), 

M jj 



(5) 



h — d 

m=l h " m 



(6) 



with the boundary conditions 



lim &(Jfe) ~ k 2l+1 and 

fe^O 



lim Se(k) ~ k . 

k— >oo 



(7) 



In any case there are 2-4 poles d m and strengths D m sufficient to provide a smooth descrip- 
tion of data. Using a [4/4] or [6/6] Pade approximation for the exponential function e 2 and 
substituting the rational phase function (pf) into z = 2i5e(k) gives a rational S-matrix 



2N " N k + al k + a}, 



Sr 

r[k - a n *=ik - ai k - an 



^) = i + ErV = nFf'Fi, (8) 



using the notation {a n } := {cr n |Im(a n ) > 0} and {cr n } : = {<7 n |Im(<7 n ) < 0}. 
The Marchenko input kernel 



1 



-oo 



F e (r,t) = -— / hj(kr) [S e (k) - 1] hj(kt)dk (9) 

Z7T 



oo 



is readily computed with the Riccati-Hankel functions hf(x) and contour integration. This 
implies an algebraic equation for the translation kernel At(r,t) of the Marchenko equation 

/•oo 

Ai(r,t) + F t (r,t) + A i (r,s)F i (s,t)ds = 0. (10) 

Jr 

The potential is obtained from the translation kernel derivative 

V e (r) = -2-^-Mr,r). (11) 
The rational representation of the scattering data leads to an algebraic form of the potential 



B. Gelfand Levitan Inversion 

Gelfand-Levitan inversion uses Jost functions as input. The latter is related to the 
S-matrix by 

Using the representation (|8|), the Jost function in rational representation is given by 

N b — n i N R 

W = n^f = i + Erx L T. (i3) 

n=l ™ i &n n=l ™ < &n 

or 
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wmr 2 = i + x:^: 

n= l ft, — a rt 



12- 



The input kernel 



2 r°° 
G e (r,t) = - / 

7T JO 



- 1 



jt(kt)dk, 



ji(x) the Riccati-Bessel functions, is analytic. The Gelfand-Levitan equation 



K e (r, t) + Gi(r, t) + f K £ (r, s)G e (s, t)ds = 0, 
Jo 



relates input and translation kernels and the potential is defined by 

V e (r)=2^-K e (r,r). 
ar 



(14) 



(15) 



(16) 



(17) 



For this potential an algebraic form is known 0. 

Gelfand-Levitan and Marchenko inversions yield the same potential. Numerical insta- 
bilities can make the potentials differ but this signals in practice a problem and thus is 
permanently checked. 



III. INVERSION POTENTIAL RESULTS 

Today, partial wave phase shift analyses are available for many hadronic systems. The 
best known of this sample is the NN analysis of Arndt et al. (SAID) which covers an energy 
range 0-1.3 GeV for np and 0-1.6 GeV for pp scattering [|1(J. The pp data are presently 
extended up to 3 GeV with measurements from the EDDA Collaboration of COSY in Jiilich 
||11|| . Of a similar quality is the nN analyses which is also available from SAID. Partial wave 
phase shifts of the nN system are determined by an analysis of elastic n + p — > n + p, ir~p — > 
n~p and charge exchange n~p — > 7r°n scattering. We use the solution SM95 of Arndt et 
al. ||12|| . This and all other analyses suppress Coulomb effects and assume good isospin. 
This implies that mass differences between and n° as well as proton and neutron are 
neglected. As a rule, for the pion is used = m^i = 139.5676 MeV and for the nucleon 
tun = m p = 938.27231 MeV respectively. Included is also the Karlsruhe-Helsinki analysis 
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of Koch and Pietarinen KH80 |13|. For tttt scattering we use phase shifts from the analysis 



of Frogatt and Petersen jnj] and theoretical predictions from a meson exchange potential 
by Lohse et al. and chiral perturbation theory by Gasser and Leutwyler ||16|| . Finally, 
the Kir analysis of Estabrooks et al. uses final state interactions of K ± p — > K ± 7T + n and 
K^p -> K ± 7T-A ++ |T7[. We restrict our analyses herein to L = and 1 partial waves 
with isospins T = 1/2 and 3/2 for ttN and Ktt, and T = 0,1 and 2 for the tttt systems. A 
comprehensive analysis of data and inversion potentials can be found elsewhere |§. We are 
using the elastic domain phase shifts and thus limit the input to < 500 MeV for ttN, 
M nrr < 970 MeV and Mk-k < 1-3 GeV. Resonance effects, like the /o(975) in tttt scattering, 
are not included. 



A. ttN s wave scattering 

The notation of this channel distinguishes Su and S31 partial waves to signal angular 
momentum L = 0, isospins T = 1/2 and 3/2 and spin S = 1/2 states. Input phase shifts 
from SM95 and KH80 are shown in Fig. 1 and the inversion potentials are given in Fig. 2 
respectively. In Sec. IV. A these potentials are used to study pionic hydrogen. 

As a brief note we mention the assessment of the pion nucleon coupling constant from 
these potentials. From the ttN potentials in the T = 1/2, 3/2 s-channels we find the 
scattering lengths a\ = 0.178 m" 1 and a 3 = —0.088 m" 1 . For a comparison with several 
other predictions see Table I. These results may be used in the Goldberger-Miyazawa- 
Oehme sum rule [EDI 



fl NN _ (ml - f){m N + m,) (fli _ ^ _ <~f ^ ^=^L dq (i 8) 



Att QuinTU^ 8tt z Jo J q 2 + m 2 

to obtain a model independent estimate for the ttNN coupling constant. Using the simplified 
form of this sum rule 



tNN 



Ait 

where the integral 



0.19m 7r (ai - a 3 ) - (0.025 mh^)J, (19) 



J 



00 a. 




=£dq 

2 



(20) 



4tt 2 Jo 



has the VPI value J = -1.041 mb. We find 




0.0766, 



or 



9nNN 
47T 



13.84, 



(21) 



which is fully consistent with the value of 13.75 ± 0.15 given in Jl2| . 



B. ttN p wave scattering 



Most prominent is the A (1232) or P 33 resonance which we treat with great care. The low 
energy phase shift function, shown in Fig. 3, uses SM95 and KH80 data. A factorization 
of the S-matrix into a resonant and a non-resonant background part S(k) = S r (k)Sb(k) is 
useful. For the resonant part a resonance and an auxiliary pole parameterization 



is used. It contains the right amount of zeros and poles for a decomposition into Jost 
functions. The background S&(fc) is smooth but its shape depends on the parameters in 



inversion potential is independent from this splitting and is shown in Fig. 4. The resonance 
feature is generated from the short range attraction near the origin and is limited by a 
narrow (0.1 fm) potential barrier with a height of ~ 20 GeV. Being accustomed to strengths 
and ranges of NN potentials it is surprising to see this small radial dimension and large 
potential strength. After a second thought this is not a surprise since inversion is a kind of 
generalized Fourier transformation and thus the units are a consequence of the pion mass 
and the quantitative behavior of the phase shift. Nevertheless, the potential changes occur 
at very small radii compared to the size of the charge form factors of pion and nucleon. 
The RMS radii are known to be 0.54 fm for the pion and 0.7 fm for the nucleon. Our 
radius describes the distance of the center of masses and r = 0.16 fm, the barrier radius, 
implying more than 90% overlap of the intrinsic structures. We conjecture for the barrier a 



S r (k) 



(k + k r )(k - k* r ) (fc + k h )(k - fcg 
(k - k r )(k + k*) (k-k h )(k + k* h ) 



(22) 



(f2"2]). Actual values are taken from data tables [f2T|] for k r and k^ = 0.5 + HO fm 1 . The 



9 



simulation of a transition from the pion-nucleon quark content into the 3-quark content of 
the A. Ultimately, such explanation must be confirmed by QCD calculations. 

In Fig. 5 we illustrate the P33 potential resonance showing the physical solution 
\u(E, r)/r| 2 as a function of energy and radius. This figure shows that the probability 
builds up between the origin and the barrier. 

The long range part of the inversion potential, not visible in Fig. 4, behaves like a 
Yukawa tail with a strength Y = 650.0 MeVfm and an effectively exchanged mass of 350 
MeV. There exists no physical particle with this mass since in the meson exchange picture 
s- and t-channel graphs contribute. 



C. 7T7T s wave scattering 

mi phase shifts come from the analysis of final state interactions in ttN — > 7T7riV systems 
or the i^-decay K~ — > 7r + 7r~ez/. Here we use results of the CERN-Munich experiment 



n|. The scattering is purely elastic until = 987.3 MeV where coupling to the KK 
opens and the phase shift becomes highly inelastic and resonant in the T = channel. The 
T = 2 channel remains smooth. A summary of all experimental phase shifts and theoretical 
predictions, from a meson exchange model [15] and chiral perturbation theory , is shown 



in Fig. 6. Notation for the isospin and angular momentum channels uses Sjor V^.The 
three used phase shift sources yield slightly different inversion potentials, shown in Fig. 
7. It is obvious that tttt phase shifts are less well established than 71N data. Table II 
comprises a summary of effective range parameters from different sources to substantiate 
the uncertainties. 



D. 7T7T p wave scattering 

As before, phase shifts come from the analysis of final state interactions in ttN — > titiN 
systems. Isospin is limited to the single value T = 1 and the CERN-Munich analysis is used 
|i"4| . The scattering is dominated by the p resonance, m p = 770 MeV, T = 150 MeV, and 
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the phase shift remains essentially real until M n7r = 1.2 GeV. The used phase shift is shown 
in Fig. 8 and the inversion potential is shown in Fig. 9. Notice the barrier maximum at 
0.16 fm which is the same as seen in the ttN P33 channel. We decline from repeating the 
resonance wave function display due to its similarity with results shown in Fig. 5. 

E. Kit s wave scattering 

Phase shifts are taken from the Estabrooks et al. analysis [fT7|| . We distinguish between 
isospins T = 1/2 and 3/2. As shown in Figs. 10 and 11, this systems resembles the nn 
s-wave scattering. The long range part of the two isospin potentials are numerically very 
close outside r = 0.4 fm which supports isospin independence. Furthermore, outside r = 0.8 
fm both potentials are very small, |V^ 2T (r)| < 1 MeV. This weak medium and long range 
interaction requires further investigations but we notice the same feature for K + N inversion 
potentials, consistent with the extraordinary long mean free path of kaons in nuclear matter 
p4],|9|. We shall show that this system displays similar medium effects as the corresponding 
7T7T system. 

IV. COUPLED SYSTEMS 

P ionic-hydrogen A np is treated in a coupled channel calculation as a charge exchange 
resonance seen in 7r°-neutron scattering. There exists a spectrum of excited states with 
different relative angular momenta but only the s-states shall be studied for comparison 
with experiment. Similarly, A nn is formed by 7r~7r + and it is treated as a charge exchange 
resonance state seen in the elastic 7r°7r° channel. 

The coupled system is written in the form 

tf + **Vi = E^Vfc/i. (23) 

3=1 n 

and uses conventionally i — 1(2) as the reaction (elastic entrance) channel. 
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A. Pionic hydrogen 



The reduced masses enter in our calculation with two options. In the first case we use 
only the mass of the charged pion and the proton mass, /ii = /x 2 with 

(m^- + m p ) 

In the second case we use /xi ^ /12 and use for p, 2 the physical masses of the neutral particles 

M2 = 7 : r- 25 

(77V + m n) 

The potential matrix contains the Clebsch-Gordan coefficients for the rotation of the inter- 
action from good isospin into particle states. 

Vu = V' = + ±V}1* - (26) 

V„ = V a = ^(VW-V}"), (27) 
V 22 = V*° n = \v?/* + \V} /2 - (28) 

The kinematics is expressed in any case by the physical masses and projectile kinetic energy 

Tlab 

S = (m n o + m n ) 2 + 2Ti ab m n (29) 
,2 _ S 2 + K_ - mlf - 2S{m%. + m}) 

,2 _ m 2 n T lab (Ti ab + 2T7V) 

k2 ~ Sihcf ' [6) 

which guarantees the correct threshold behavior. 

The Coulomb attraction between ir~ and p causes a bound system which is known as 
pionic hydrogen A np . In addition to Coulomb attraction, the hadronic interaction between 
the two constituents distorts the short range interaction and changes the pure Coulomb 
spectrum. Decay channels are ir~p — > 7r°n + 3.30023 MeV and n~p — > 727. The hadronic 
shift of the 3p — > Is transition and the total Is width has been measured at PSI ||. To 
analyze this experiment we use the described potential model. Inelasticities, Coulomb and 
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other isospin breaking effects are supposed to be not included in the SM95 phase shift 
analysis and the real potential matrix is shown in Fig. 12. We use the reduced masses 
l^i — A*2 = 121.4970 MeV which are consistent with the phase shift analysis and using 7r ± 
and p masses. Alternatively n\ = 121.4970 MeV and fi2 = 118.0216 MeV based upon 7r° 
and neutron masses respectively can be used. 

The bound states of the ir~p system can be found as resonances in the energetically 
open 7r°n channel. The width (FWHM) of this resonance accounts only for the decay of 
7i~p — > 7T°n. Table III contains this result in the first line. To evaluate the hadronic shift 
we used in all cases the reference energy Eg = 3234.9408 eV. The shift is calculated as the 
difference between Eg and the calculated resonance energy taken at the maximum of the 
elastic scattering cross section 

a{n°n^Tr Q n) = ^\l-S Xi \ 2 . (32) 
fc 2 

From this distribution we obtain also the FWHM. The Coulomb potential is contained in V%\ 
as determined from point charges V = e 2 /r or double folded Gaussian charge distributions 
V° = e 2 $(1.13r)/r. $(r/a) is the error function and a = ^J< r 2 > + < r 2 > with RMS 
radii < r 2 > 1 / 2 = 0.5389 fm and < r 2 > l ' 2 = 0.702 fm. To compare with experimental data 
H we evaluated the partial width T'fj 1 using the Panofsky ratio P = 1.546 ± 0.009 |35| , 

T=(l + ^j T^ p ^° n . (33) 

To account also for decay into the channel from the 7r°n channel an imaginary potential 
W\\ = — 9exp(— 9r 2 ) MeV is added to Vn which brings the theoretical results in agreement 
with experiment. Other corrections are not further pursued since they enter into the refer- 
ence and resonance energy with approximately the same amount, affecting the difference by 
< ±0.01 eV. We shall not dwell upon this issue in more detail herein. The essential result of 
this part of our calculation is that isospin breaking effects due to mass differences are small. 
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B. Pionium 



Very close the same convention as for A np is used here. Only one value of a reduced mass 
occurs 

= 1*2 = (34) 

This value is varied in Sec. V.A when medium effects are discussed. The potential matrix 
is readily expressed by 

Vn = = \V> + 2 -V° - ^ (35) 



V 12 = V 21 = ^(V 2 -V °), (36) 

2 „ 1 
3* + 3 



V* = V' v = %Vl + V (37) 



and the kinematics by 



S = Amlo + 2T lab m^ (38) 

, 2 _ ml Ti ab (Ti ab + 2m^>) 
M m = y/S (41) 

Similar to the A np system there exists pionium A n7T which is formed by 7t~tt + Coulomb 
attraction. It decays predominantly by charge exchange into the open 7r°7r° channel. The 
coupled channel system is defined with Eqn. (p3|). We assume the same approach as for A np 
and rotate the good isospin potentials into particle states. The hadronic potential matrix 
is shown in Fig. 13 using the three different sources discussed in Sect. III.C. Phase shift 
analyses and inversion use a single mass ^1 = ^2=/^ = m n +/2 without Coulomb effects. 
This assumption guarantees good isospins T = and 2. The results of our calculations are 
summarized in Table IV. A point Coulomb reference energy E± s = 1.85807248 keV is used. 
With this choice of masses and the correct Q-value our calculations agree practically with 
some well known scattering length expressions which is most often used |6|, j36| , [37H 
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l_87r/2Amy (gjj - |^(0)| 2 , . 

r 9 ^ // J l + |/iAm(a[l + 2a2) 2 

« 1.43 (a°-a 2 ) 2 |vD( )| 2 , (43) 



where Am is the mass difference m n + — m n o . The real interesting result of this study comes 
from dramatic changes (shortening) of the A nn lifetime, by several orders of magnitude, with 
small variations of the reduced masses /ij. 

V. MEDIUM EFFECTS 

Low energy nuclear physics associates medium effects with changes of the free interaction 
potential or free two body scattering amplitudes in the presence of a few- or many-body 
environment. Effective interactions incorporate such effects and degrades the environment 
into the role of spectators. A well known example is the nuclear matter g-matrix which 
includes Pauli blocking and selfconsistent mean field effects and which is related to the free 
two body t-matrix. Other medium modification are due to recoil, truncation of coupled 
channels and relativistic corrections to name a few. Medium effects from meson and boson 
exchange models are related to restoration mechanism of the broken chiral symmetry in 
nuclear matter. This field of research is presently in the center of different theoretical lines 
of thought and some of them are found in [[38] . 



Our concern is an effective mass in the two particle wave equation which may have dif- 
ferent causes in few and many body systems. At this stage the consequences and dynamical 
effects are shown for the two particle subsystems tttt and Ktt in case an effective mass is 
assumed. What we study are partial wave phase shifts S(E,m) as a function of energy and 
the effective mass, which determines the reduced masses //j in Eqn. (0). Effects upon the 
lifetime of pionium are also studied. 
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A. 7T7T scattering 



In Sect. IV. B results of pionium lifetime and hadronic shift are given. This study uses 
the tt^ mass consistently with the phase shift analyses and the scattering length expressions 
for the life time. Here we extend this study and show results in Fig. 14 for the hadronic shift 
and lifetime of pionium as a function of the pion mass, m n = 2/xi = 2/i 2 in Eqs. fl3~4"D. The 
potential matrix Vij was unchanged and the kf values computed with the physical masses 
of 7T and 7T° as defined in Eqs. (38-41). With this prescription we effectively change the 
total strength of the potential matrix by ±1% and observe a variation of the shift by two 
orders of magnitude and of the lifetime by three orders of magnitude. The physical mass 
result is accentuated. The interval of mass variation is small compared with the physical 
mass differences between 7r ± and 7r° of 4.6 MeV, which a correct phase shift analysis should 
consider. Since this is not the case, consistency requires to use fii = fi2 in the inversion 
procedure and in the coupled channel calculations. Thus, we maintain the assumption of a 
good isospin also for the calculation with an effective mass. 

Triggered from this result we extended our calculations for the eigenchannel phase shifts 
of the coupled system which coincide with the uncoupled good isospin calculations for T = 
and 2. In Fig. 15 we show the L = 0, T = phase shift 5q(to, T) as a function of the 
effective mass m = 1\x\ = 2/z 2 and the laboratory kinetic energy. The physical mass result 
is emphasized and follows the result shown in Fig. 6 for 5$. As a function of mass we 
observe a typical resonance behavior whose width T = 1 — 2 MeV which is almost energy 
independent. The resonance feature is supported by Fig. 16 which shows \u(r, m, T)/r\ 2 as 
a function of radius and effective mass for three kinetic energies. Since L = 0, the radial 
wave function is 7^ at the origin around which the probability is build up. This figure 
should be compared with Fig. 5 of the nN A resonance which shows the same pattern as 
a function of energy. For the ttN P33 phase shift 5(E) varies only insignificantly when the 
reduced mass of pion-nucleon system is modified within 1 MeV and the resonance pole is 
stable against this variation. From the nN potential in Fig. 4 the same mass dependence 
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as for the irn system was expected. This is not the case. We attribute this stability to the 
boundary condition of the p-wave at the origin to be zero. The tttt s-wave T = potential 
in Fig. 7 supports the resonance as function of mass since no boundary condition at the 
origin restricts the development of a large amplitude breathing mode. The potentials have 
in both cases a very high barrier of 20 - 30 GeV, compared to the kinetic energy of several 
100 MeV, which is narrow ~ 0.1 fm and within the barrier the potentials are very deep. For 
the s- and p-waves decays the radial wave exponentially practically independent from the 
projectile energy. However, the boundary condition of the p-wave at the origin suppresses 
the free unfolding of a resonance amplitude independent from energy. It requires an optimal 
matching of the external wave function to realize the internal resonance enhancement. A 
small variation of the potential depth within the barrier does not overcome the restriction 
from the boundary condition for the p-wave. This situation is essentially different for s- 
waves. For them it is important to have the correct wave number within the barrier that the 
wave function matches the exponentially decaying function in the barrier optimally. The 
wave number within the barrier is determined by 

A small variation of the reduced mass n produces the discussed effect. In connection with the 
very high barrier gives this the explanation of the projectile independence of the observed 
7T7T resonance as a function of the effective mass. The high barrier decouples for s-waves the 
inner from the exterior dynamics which is not the case for p- and higher partial waves. 

Finally, the effective mass used in this study must be the result of embedding the nn 
system into a few- and many-body environment. In the meson exchange models for NN 
scattering should the correlated two pion exchange dynamics not be determined from the 
physical free pion mass but from an effective mass. This mechanism is able to change the 
non-resonant two pion system into a resonant system with a width 1 < T < 600 MeV. The 
lower limit is extracted from Fig. 15 and the upper limit is deduced from the free system 
<5q (e) |§. This gives a possible explanation of the properties of the a-meson used in all high 
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quality meson echange model for hadron-hadron interactions out of which the NN potentials 
are the best examples. The phase shift 8q{E,w) shows practically no mass dependence and 
maintains its repulsive nature. 

To support this picture for scattering it is obvious to look for other systems with the 
same properties. 

B. K-k scattering 

Phase shift data for the K7r s-wave scattering yield qualitatively the same situation with 
two possible isospin couplings T = 1/2 and 3/2. The inversion potentials are shown in Fig. 
11. A purely repulsive potential is seen for the T = 3/2 channel whereas T = 1/2 has a 
narrow and high potential barrier at about the same radial region as the tttt system and we 
expect great similarities with respect to the phase shifts 5g(m, T), a functions of effective 
mass and kinetic energy. Fig. 17 confirms this conjecture. Following the same line for 
the coupled system in Sec. IV, one effective mass parameter is used m = I^k-k where the 
physical mass has a value of 215.94 MeV. The solid line signals the physical T = 1/2 phase 
shift and the energy-mass distribution shows a resonance structure. The width increases, 
T > 80 MeV, with laboratory kinetic energy which is caused by the relative smallness of the 
barrier of 800 MeV. The phase shifts 5q(E, m) shows no mass dependence and maintains its 
repulsive nature. 

VI. SUMMARY 

The elastic scattering domain for ttti, Kit and ttN scattering is investigated with the 
help of a local r-space potential model. Quantum inversion is used to generate from phase 
shift data the potentials for low partial waves and the permitted isospin channels. Tables of 
these potentials are available |39j . All resonance features visible in the phase shifts are well 
accounted for, p and A are potential resonances of a tctt and tcN system respectively. 
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Pionic hydrogen A np is studied as a resonance in elastic 7r°-neutron scattering below 
7r~-proton threshold. The calculated shift and width of the ground state agree very well 
with recent measurements. Similarly, pionium A n7T is studied as a resonance in elastic 7r°7r° 
scattering below 7r + 7r~ threshold. The potential model confirms well known scattering length 
expressions often used with hadronic bound state problems. As a new feature we disclose 
a medium effect of lifetime and shift of A n7T when the interaction strength is increased by 
typically half a percent. The amount of change is more than three orders of magnitude. In 
an extended study a medium resonance is identified with the a resonance. In particular it 
is shown that this resonance is a function of the effective mass, which we identify as the 
factor in front of the potential, and not as it is common as a function of energy. It can 
be considered as a parametric resonance. The change of A n7r lifetime is a feature of this 
resonance. A similar situation with medium effects and a parametric resonance is identified 
in the Kix L = 0, T =1/2 channel. Contrary to this medium resonances show the classic p 
and A resonances practically no effective mass dependence. Showing radial wave functions 
we explain these differences of dynamical behavior with boundary conditions at the origin. 
All resonances are associated with a small radial region within 0.3 fm. This distance is 
the center of mass distance of the two particles and is not the QCD bag radius. Table V 
shows that the potential barriers which trap the relative system are very close the same 
with the exception of the strange Kir resonance. In view of the large RMS radii of mesons 
and nucleons (~ 0.6 fm) this implies a more than 90 percent overlap of the individual QCD 
bags before fusion occurs and the new configuration develops. We see in this universal small 
barrier radius a reason why meson exchange works and permits quantitatively to describe 
NN and other hadron-hadron interactions. This finding requires detailed calculations on 
the QCD level. 

For low and medium energy nuclear physics this potential model defines essential ingredi- 
ents of boson exchange models in terms of 7r mesons and nucleons only, p and A resonances 
are well accounted for their properties with the inversion potentials and the a meson shows 
dynamic aspects which does not support it as a particle. Since long there exist experimental 
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data which involve two pion production and still lack of satisfactory explanation. It is known 



as ABC effect |^D| and we suggest a re-analysis of this data with the discussed dynamical 
effects included. The set of inversion potentials can easily be completed || to have a basis 
for potential model calculations which use only pions, protons and neutrons. This opens the 
possibility for quantitative fragmentation calculations requiring massive supercomputing. 
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TABLES 

TABLE I. ttN s-wave scattering lengths and ttNN coupling constant obtained from the GMO 
sum rule. 



Model 



03 ] 



Ref. 



SM95 Inversion 
KH80 
ir~p Is state 
Pearce et al. 
Schiitz et al. 



0.178 
0.173 
0.185 
0.151 
0.169 



-0.088 
-0.101 
-0.104 
-0.092 
-0.085 



0.0766 
0.079 
0.081 
0.072 
0.074 



HI 
I 

El 




TABLE II. 7T7T s-wave scattering lengths. 



Source/Model 



«0 Kr 1 ] 



K - °o) K 



Ref. 



Weinberg 



Predictions from theory 
0.16 -0.046 



0.206 



XPT 


0.20 


-0.042 


0.242 


Meson exchange 


0.31 


-0.027 


0.337 




Results from experiment 




K e4 


0.26 ± 0.05 


0.028 ± 0.012 


0.288 ±0.051 


Chew-Low PSA 


0.24 ±0.03 


-0.04 ± 0.04 


0.280 ± 0.050 


Soft-Pion 


0.188 ±0.016 


0.037 ± 0.006 


0.225 ±0.017 


ttN — > ttttN 


0.205 ±0.025 


0.031 ±0.007 


0.236 ± 0.026 




Results from 


inversion 




Froggatt phase shifts 


0.31 


-0.059 


0.369 


XPT phase shifts 


0.20 


-0.043 


0.243 


Meson ex. phase shifts 


0.30 


-0.025 


0.325 



n 



El 

m 
m 
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TABLE III. A wp Is level shift with respect to Eg = 3234.9408 eV and width. The strength of 
the imaginary W\\ was adjusted to reproduce the experimental value. 





Point Charge Coulomb 


Gaussian 


Charge Coulomb 


Mass 


Shift [eV] FWHM [eV] 


Shift [eV] 


FWHM [eV] 



Mi = M2 -7.29821 0.5250 -7.16746 0.5230 

Mi + M2 -7.13259 0.5187 -7.01055 0.5144 

Shift = -7.127±0.046 [eV], Yf s n = 0.590 [eV] Sigg @ 

Mi + M2 -7.23259 0.9763 -7.12387 0.9763 V n + iW n 

Shift = -7.127±0.028±0.036 [eV], T ls = 0.97±0.10±0.05 [eV] Sigg [| 



E ls [keV] 



TABLE IV. A n7T properties from inversion potentials 



Shift [eV] 



r [10 



-15 



FWHM [eV] 



Ref. 



1.8638814 
1.8635114 
1.8616174 

1.858 
1.865 



-5.809 
-5.439 
-3.545 



1.97 
1.89 
3.22 



0.3481 
0.3627 
0.2128 



Predictions from experimental analysis and other models 



-7.0 



2.91& 



3.2 



Froggatt |L4| 
Lohse 111 

XPT 

Afanasyev | 
Efimov |3l 



TABLE V. Potential barrier positions. 



Channel 


vriV-A 


mr-a 


TTTT-p 


Ktt-K$ 


Barrier radius [fm] 


0.165 


0.145 


0.152 


0.303 



25 



FIGURES 
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3 r 




FIG. 5. ttN P33 channel radial probability distribution as function of energy, Ti^. 
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1.0 



FIG. 6. tttt L = 0, T = and 2 phase shifts. Froggatt data [14] and interpolation (dots and 
solid line), xPT (l(| (crosses), Estabrooks et al. p^j (boxes), Grayer et al. (diamonds), Manner 
(triangles) and Baillon et al. pll (asterixes). 
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0.0 0.5 1.0 1.5 

M [GeV] 

7T7T L J 

FIG. 8. 7T7T L = 1, T = 1 phase shifts. Data from Froggatt et al. @ (dots), Ochs (squares), 
Protopopescu et al. |[|] (triangles) and Deo et al. |3j| (asterixes). Inversion result (solid line). 
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FIG. 15. 7T7T L = 0, T = phase shift as a function of energy (T) a &) and effective mass. The 
physical value m n + = 139.5676 MeV is emphasized. 
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100 MeV 




FIG. 16. 7T7T L = 0, T = radial probability distribution as a function of effective mass given 
for three energies {T lab = 100, 500 and 900 MeV). The physical value = 139.5676 MeV is 
emphasized. 
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FIG. 17. Kit L = 0, T = 1/2 phase shift as a function of energy (T[ a b) and effective mass. The 
physical value 2/j.Kw = 215.94 MeV is emphasized. 
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